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ABSTRACT 

We study the nature of faint, red-selected galaxies at z ~ 2 — 3 using the Hubble Ultra Deep 
Field (HUDF) and Spitzer Infrared Array Camera (IRAC) photometry. Given the magnitude limit 
of the HST data, we detect candidate galaxies to Hab < 26 mag, probing lower-luminosity (lower 
mass) galaxies at these redshifts. We identify 32 galaxies satisfying the (Jno — i^i6o)AB > 1.0 mag 
color selection, 16 of which have unblended [3.6/im] and [4.5/im] photometry from Spitzer. Using this 
multiwavelength dataset we derive photometric redshifts, masses, and stellar population parameters 
for these objects. We find that the selected objects span a diverse range of properties over a large 
range of redshifts, 1 ^ ^ ^ 3.5. A substantial fraction (11/32) of the (Jno — i^i6o)AB > 1 mag 
population appear to be lower-redshift {z < 2.5), heavily obscured dusty galaxies or edge-on spiral 
galaxies, while others (12/32) appear to be galaxies at 2 < z < 3.5 whose light at rest-frame optical 
wavelengths is dominated by evolved stellar populations. We argue that longer wavelength data 
(> 1 /im, rest-frame) are essential for interpreting the properties of the stellar populations in red- 
selected galaxies at these redshifts. Interestingly, by including Spitzer data many candidates for 
galaxies dominated by evolved stellar populations are rejected, and for only a subset of the sample 
(6/16) do the data favor this interpretation. These objects have a surface density of ~ 1 arcmin"^. 
We place an upper limit on the space density of candidate massive evolved galaxies with 2.5 < z < 3.5 
and i^i6o(AB) < 26 mag of n = 6.6I3 Q x 10~^ Mpc~^ with a corresponding upper limit on the stellar 

mass density of p* = 5.6I2 8 ^ ^0 Mpc~^. The z > 2.5 objects that are dominated by evolved 
stellar populations have a space density at most one-third that of z ~ red, early- type galaxies. 
Therefore, at least two-thirds of present-day early-type galaxies assemble or evolve into their current 
configuration at redshifts below 2.5. We find a dearth of candidates for low-mass (< 2 x 10^^ M©) 
galaxies at 1.5 < z < 3 that are dominated by passively evolving stellar populations even though the 
data should be sensitive to them; thus, at these redshifts, galaxies whose light is dominated by evolved 
stellar populations are restricted to only those galaxies that have assembled high stellar mass. 
Subject headings: cosmology: observations - galaxies: evolution - galaxies: high-redshift - galaxies: 
stellar content - infrared: galaxies 



1. INTRODUCTION 

In the hierarchical picture of galaxy evolution, galaxies 
asse mble within dark matter haloes that merge oy er time 
(e.g. , iFreedman eT^ l2Q0ll : ISpergel et all 120031 ) . This 
model has been highly successful at reproducing galaxy 
clustering properties based on well-understood gravita- 
tional physics (e.g., Weinberg et al. 2004). However, it 
is difficult to test how galaxies assemble their stellar con- 
tent within the context of these models because of un- 
certain physics coupling the dark matter and baryonic 
matter in these haloes. 

Observationally, the cosmic star-formation-rate (SFR) 
density rises by roughly an order of magnitude from 
z ^ to 1 and remains roughly constant from z ^ 1- 
6 (|Hopkinsl l2004l . and references therein). At the same 
time, roughly half of the stellar-mass density in galax- 
ies formed between z ~3 and 1, corresponding to this 

^ This work is based in part on observations taken with the 
NASA/ESA Hubble Space Telescope, which is operated by the As- 
sociation of Universities for Research in Astronomy, Inc. (AURA) 
under NASA contract NAS5-26555; and on observations made with 
the Spitzer Space Telescope, which is operated by the Jet Propul- 
sion laboratory, California Institute of Technology, under NASA 
contract 1407. 

^ Spitzer Fellow 



peak in the SFR d e nsitv fe.^.. [Dickinson et al.l l2003l : 
lR,udnick et all [20031 : iGlazebrook [200I . While 

theoretical models broadly reproduc e this evolution 
in the global galaxy population (e.g., ISomerville et al.[ 
I2OOII : iHernquist fc Springell [2003), they have diffi- 
culty in reproducing star-formation trends in individ- 
ual galaxies, particularly in accounting for the forma- 
tion and evo lution of specific SFRs of the most mas- 
sive objects (iBaugh et al. 2003; Somer ville et all l2004l : 
[De Lucia et al.[[2006[ : [Croton et al. 2006). In particular, 
it is difficult for models to produce galaxies (at any red- 
shift) that are mostly devoid of current star-formation 
and e volve only passively (e.g.. ^ Croton et aT]l2006l : fPavel 
[2005[ ). Quantifying the number of passively evolving 
galaxies — and as a function of redshift and stellar mass 
— is helpful for constraining the processes that form stars 
within galaxies. 

Searches for high-redshift galaxies that are dominated 
by older, passively evolving stellar populations are chal- 
lenging. The light from these stars peaks at rest-frame 
optical and near-IR wavelengths (:^ 0.4 — 2 /im), which 
at z > 1 are shifted into the observed-frame near-IR. 
Several surveys have used deep, near-IR observations to 
search for galaxies up to z ~ 1 — 2 whose light_is_dom- 
inated by older stellar populations (see e.g., IMcCarthv[ 
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l2QQ4f ). Such objects should have very red R — K colors, 
which span the 4000A/Balmer-break at these redshifts; 
several searches have focused on so-called extremely red 
objects (EROs), with selected objects typically satisfy- 
ing {R — K)yQgg, > 5 mag. As a class, EROs include 
both pas sively ev olving early-type galaxies up to z ~ 1 
(e.g., Ci matti et a l. 2002) and dust-reddened starbursts 
(e.g.. Small et al. 2002). However, with only optical and 
near-IR photometry, it was not possible to uniquel y inter- 
pret t he majority of thes e objects (e.g., M oustakas et al.l 
I2004J : IWilson et aT]l2007l ). Recent observations at longer 
wavelengths from Spitzer indicate that more than half 
of the E RO popula tion have strong emission in the ther- 
mal IR (jWilson et al. 2004; Yan et al. 2004b), suggesting 
that such objects constitute a significant portion of the 
red-selected galaxy population. 

Recently, observers have used deep near-IR data 
to study r ed-selected galaxies at higher redshifts. 
iFranx et all (2003) used deep VLT/ISAAC JHKs ob- 
servations to identify a population of "distant red galax- 
(DRGs) with {Js - Ks) Vega > 2.3 mag, which 



les 



should select galaxies that have a strong Balmer/4000 A 
break betwe en the J . and Kc, bands at z ~ 2-3.5 



(see also Saracco et a l.l 120011). Like EROs, the DRG 
color selection is also sensitive to high redshift {z > 
1) galaxies dominated by dust-enshroude d starbursts 
also have red observed J — K colors (e.g., ISmail et al.l 
l2002f ). Indeed, subsequent study of DRGs has shown 
that the majority are massive gal axies that are ac- 

1.5-3 (Ivan Dokkum et all 



tive_ly forming stars at z 

1 iForster Schreiber et all l2004 
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R.ubin et all 2004; 



Knudsen et al . 2005; Red dv et al.]l2005l : iPapovich et"al] 
20061 : IWebb et al.ii2006) . and only a small subset appear 



to be con ipletely devoid of star formation and passively 
evolving (iLabbe et al1l2005l : iKriek et al.ll2006f ). 

However, nearly all surveys of red-selected objects have 
been restricted to using near-IR images from ground- 
based telescopes, which in practice limits these searches 
to i^AB ^ 24 mag (e.g., Labbe et al. 2003). For passively 
evolving stellar populations at z ~ 3, this magnitude 
limit acts as a limit in stellar mass of > 10^^ Mq. In 
contrast, the z ~ 2 — 3.5 UV-selected, star-forming galax- 
ies (so-called Lyman-break g alaxies. LBG s) that domi- 
nate the SFR density (Reddv et al.l [20051 ) have typical 
stellar masses - 10^^ Mq (for "L*" LBGs; Papovich et 
al. 2001). T he inferred ages for LBGs a re generally less 
than 1 Gyr (|Shaplev et all l200lL l2005f ) — significantly 
less than the age of the Universe at these epochs. Thus, it 
is unknown if there exists a population of "faded" LBGs, 
identifiable as red, passively evolving galaxies with stellar 
masses ~ 10^^ Mq. Deriving constraints on the density 
of such objects would improve our understanding of the 
evolution of these galaxies and on how galaxies assem- 
ble their stellar content at these early epochs. However, 
detecting lower-mass, red galaxies at high redshift re- 
quires near-IR surveys with higher sensitivity than what 
is practically available from the ground. 

Here, we use observations in the H ubble Ultra Deep 
Field (HUDF; Beckwith et al. 2006 and lThompson et al.l 
12005). This field has the deepest optical HST/ ACS im- 
ages to date, which combined with extremely deep, near- 
IR HST /memos, and Spitzer (3.6, 4.5 /im) datasets 
form part of the Great Observatories Origins Deep Sur- 



vey (GOODS). The depth of the images in this field pro- 
vides a means to explore high-redshift galaxies to a rel- 
atively low limiting stellar-mass. Here, we focus on a 
sample of galaxies selected with red Jno — i^ieo colors 
in an effort to identify candidates for p assively evolving 
galaxies at z > 2. Similarly, Bram mer fc van DokkunJ 



(|2007l ) use 8i J — H color selection to identify red galax- 
ies at 2: > 2 using the deep FIRES near-IR survey 
imaging. The HST /mCMOS data for the HUDF allow 
us to probe objects with lower stellar masses than has 
been p reviously possible. For example, Dickinso n et al.l 
(|2000r ) used deep NICMOS d ata in the Hubble Deep Field 
(HDF; I Williams et al. I [19961 ) to study an unusually red 
(^110 — HiQQ > 2 mag) source, possibly a passively evolv- 
ing z ^ 3 galaxy with M > 3 x 10^^ Mq. Similarly, 
Chen & Marzke (2004) used the HUDF HST images to 
identify red galaxies (^775 — i^ieo > 2) with photomet- 
ric redshifts Zph > 2.8 (see also Yan et al. 2004). In this 
work we use the HUDF data to study the number density 
and mass density of objects at 2.5 < z < 3.5 with colors 
consistent with older passively evolving stellar popula- 
tions. To study the HUDF Jno — i^ieo > 1 sources, 
when possible we include GOODS IRAC 3.6-4.5 /im and 
MIPS 24 /im data, which allows us to better constrain 
the nature of the stellar populations in these objects (as 
we discuss below). 

Throughout this paper we assume the following cos- 
mological parameters: Qm = 0.3, Qa = 0.7 and 
Ho = 70 km s~^ Mpc ~^ (consistent wit h the lat- 
est WMAP results, e.g., "Spergel et_alj [20031) . We use 
^435^606^775^850-^110^160 to denote magnitudes derived 
in the the HST ACS and NICMOS bandpasses F435W, 
F606W, F775W, F850LP, FllOW and F160W respec- 
tively. We denote magnitudes derived from the IRAC 
channel 1 and channel 2 data as [3.6/im] and [4.5/im]. For 
reference, the effective wavelengths of these bandpasses 
are 0.43, 0.60, 0.77, 0.91, 1.1, 1.6, 3.6 and 4.5 /im. Ah 
magnitudes henceforth are on the AB system, m(AB) = 
23.9- 2.5 log(/,/l/iJy). 

2. RED GALAXIES IN THE HUDF: PHOTOMETRY AND 
SELECTION 

We combine HST/ ACS and NICMOS observations 
of the HUDF (Beckwith et al. 2006^; Thompson et 
al. 2005) over a 2.4 x 2.4 arcmin^ area with Spitzer 
IRAC observations from GOODS (M. Dickinson et al., 
in preparation)^. The plate scale of the drizzled images 
is 0.03, 0.09 and 0.6 arcsec pix-^ for the ACS, NIC- 
MOS, and IRAC data respectively. The lOa AB point 
source detection limits for the ACS and NICMOS im- 
ages are 29.2, 30.0, 29.7, 28.7, 27.0 and 27.0 mag for 
the ^435 "1^606^775^850-^110^160 images respectively. For 
the IRAC bands, the 5a AB point source detection lim- 
its are 26.3 mag (0.11/iJy) and 25.6 mag (0.2/iJy), for 
[3.6/im] and [4.5/im] (M. Dickinson, private communica- 
tion, 2005). 

We use ACS images reprocessed by [Thompson et al.l 
(|2005r ) and matched to the NICMOS plate-scale and ori- 
entation. We then convolve the ACS images by a kernel 
to match the point spread function (PSF) of the Hiqq- 
band image. Because the PSFs of the Jno and i^ieo im- 
ages are similar we do not convolve the Jno-band image. 

^ http://www.stsci.edu/hst/udf 

^ http:/ /www. stsci.edu/science/goods/ 
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Fig. 1. — Color-color plot of all catalog objects with ifieo < 26.0. 
The Jii o — if 160 > 1-0 color cut is marked. The color evolution 
of three IBruzual CharlotI (|2003l ) models is shown, with integer 
red-shift steps marked. The curve labeled SSP (orange) is a l.OGyr 
old simple stellar population formed in a burst with no subsequent 
star formation. The CSP curve (blue) is a lOOMyr constant star- 
forming model attenuated using Calzetti et al. (2000) dust. The 
purple curve a is superposition of these two SSP and CSP models, 
with a ratio (SSP:CSP) of 99 : 1 by mass. Although the SSP and 
CSP models do not individually characterize the data well, the two- 
component model does. The median error for the Jno — i^ieo color 
is smaller than the symbol size; the median error for the Vqqq — zsso 
color is similar to the symbol size. 



We use standard IRAF packages to perform the image 
convolution. Our tests show that point-source photome- 
try on the convolved images recovers the same fraction of 
the total flux to within 5% for aperture radii greater than 
5 pixels (0.45 arcsec). Source detectio n was executed us- 
ingthe SExtr actor software (v2.3.2, Bertin fc Ar nout"sl 
Il996f ) on a summed Jno + i^ieo image, and photometry 
was performed on the individual PSF-matched images. 
We use the same detection and photometry parameters 
as in [Thompson et aH (|2QQ5l ). We restrict our analysis to 
objects with i^ieo < 26 mag measured in the SExtractor 
elliptical, quasi-total ("Kron"; SExtractor MAG_AUTO) 
apertures. The NICMOS data detect galaxies to 27th 
magnitude, and adopting a i^ieo < 26 mag limit en- 
sures that we can derive robust colors for objects with 
^110 - Hieo = 1. We measure ACS and NICMOS colors 
for objects in the PSF-matched images using isophotal 
(MAG-ISO) apertures derived from the i^ieo-image; we 
correct the isophotal magnitudes to total magnitudes us- 
ing the difference between the i^ieo-band MAG_AUTO 
and MAG-ISO magnitudes. 

The photometric errors from SExtractor are underes- 
timated because they do not take into account pixel-to- 
pixel correlations from, for example, the drizzling and 
PSF-convolution processes. We therefore estimate pho- 
tometric errors from the binned, drizzled, PSF-convolved 
ACS images themselves. We measure the sky noise as a 
function of aperture size in blank regions of each im- 
age. We then parameterize the noise as a power law, 
a = A X rip-^, where a is the noise measured in an 
aperture with a size of npix pixels, and we flt for P 
and A. In order to measure the noise in the convolved 
images we tabulate 1250 randomly positioned aperture 
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Fig. 2. — Color vs. redshift diagram of lBruzual CharlotI (l2003h 
population synthesis models showing Jiio— ^160 color evolution for 
different ages, star-formation histories and dust attenuations. The 
solid curves show a simple stellar population (SSP) model at two 
different ages (1.0 Gyr and 0.5 Gyr). The dotted curves show a 100 
Myr constant star-formation model (composite stellar population: 
CSP) with two different dust attenuations: E(B — V) = 0.3 (lower 
dotted curve) and E(B — V) = 0.6 (upper dotted curve). The 
dot-dashed curve (labeled SSP+CSP) shows an example of a two- 
component model: an old population (1.0 Gyr SSP) and a young 
population (100 Myr constant star-formation with E(B — V) = 0.3) 
are combined with a ratio (SSP:CSP) of 99 : 1 by mass. The dashed 
curve shows a maximally old SSP formed at z = 7. Note that our 
selection criterion, Jno — Hiqq > 1.0, generally excludes objects 
at redshifts less than z ~ 1.5. 

fluxes; we then bin these data and fit for the width of 
the distribution with a Gaussian function. We exclude 
bins in the non-Gaussian positive tail of the distribu- 
tion (these correspond to apertures that receive some 
flux from bonafide sources in the image). This proce- 
dure is repeated for 11 aperture diameters in each band: 
0.54, 0.60, 0.80, 0.99, 1.15, 1.30, 1.40, 1.53, 1.65, 1.80 and 
2V25. We obtain best-fit values of f3 = 0.74, 0.92, 0.96, 
0.93 and A = (5.6, 3.5, 2.0, 2.2) x -lO""^ ADU sec-\ for 
^435^06^775^850 respectively. Using these fitted values, 
we can then derive the error, a in each band for any 
isophotal aperture size, npix- 

In figure [H we plot the Jno — i^ieo vs. Vqqq — Zs50 col- 
ors for all HiQo < 26.0 catalog objec ts. We a lso show the 
color evolution of three Bruzual & CharlotI (2003) mod- 
els, a passively evolving burst (SSP) model, a constant 
star-forming (CSP) model, and a composite SSP+CSP 
model, at fixed age and as a function of redshift. The 
composite model characterizes the Jno — i^ieo > 1 se- 
lected objects well. 

To target galaxies whose light is dominated by evolved 
stellar populations, we select objects with Jno — i^ieo > 
1.0 mag and i^ieo ^ 26 mag, which in principle se- 
lects objects with a prominent 4000 A/Balmer break at 
z > 1.5. In figure [2] we plot se veral fiducial stella r popu - 
lation synthesis models from IBruzual fc CharlotI (|2QQ3[ ), 
with solar metallicity and a Salpeter IMF. Our color cri- 
terion excludes model stellar populations with recent star 
formation (and low-to- moderate dust attenuations), and 
should be sensitive to passively evolving, older, simple 
stellar populations (SSPs), or dust-enshrouded objects 
with younger stellar populations. We find 32 objects sat- 
isfying Jiio — HiQo > 1.0 mag and i^ieo ^ 26.0 mag in 
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Fig. 3. — i^ieo — [4.5/im] versus 2775 — ifieo colors for the 16 
galaxies in the Jiio — ^160 > 1 mag sample with robust IRAC pho- 
tometry. The lines show colors of various model stellar populations 
formed at 2; = 7. The dashed line shows the expected colors of a 
passively evolving stellar population observed at z = 2.7, formed 
in a burst as a function of age, the dot-dashed line shows the 
same model observed at z = 3.5. Tick marks denote the model 
ages along the dashed and dot-bashed curves as labeled. The solid 
lines show colors of a stellar population with constant star forma- 
tion formed at z = 7 as a function of age and extinction with color 
excess E{B — V) = 0.2 and 0.6, as labeled. The encircled data 
points indicate the 6 candidates for distant galaxies dominated by 
old stellar populations (see § 4.1). These models roughly bound 
the galaxy colors, and indicate that including the IRAC data we 
can understand the properties of the stellar populations in these 
galaxies. Note that none of the objects have colors consistent with 
passively evolving stellar populations older than 1 Gyr. 

the NICMOS-selected catalog; their properties are sum- 
marized in table [TJ 

We match Jno — i^ieo > 1 mag objects against the 
Spitzer /IUAC catalog from the Great Observatories Ori- 
gins Deep Survey (GOODS; M. Dickinson et al. in prepa- 
ration). Source detection for this catalog was performed 
on a summed [3.6/im] + [4.5/im] image, and photometry 
was measured in circular, 4.0''-diameter apertures, which 
we correct to total magnitudes by adding —0.36 and 
—0.31 respectively to the [3.6/im] and [4.5/im] magni- 
tudes (appropriate for GOODS, M. Dickinson et al., in 
prep, and slightly larger than those in the IRAC data 
handbook^). Because the image PSF size is substan- 
tially different between the NICMOS and IRAC images 
(FWHM 0'.'3 and 1'.'6, respectively), we do not attempt 
to PSF match these images. 

We matched these objects to sources in the 
Spitzer /IUAC catalog first using a 4''-radius criterion. 
We then visually inspected each source and accept only 
those objects whose positions coincide with the HST 
source (they are well-centered) and who are isolated (un- 
blended) from other HST sources within that radius. 
This is a conservative approach, and has the advantage of 
limiting our analysis to only those Spitzer counterparts 
that are free from source confusion. Of the 32 objects 
in the parent sample, 16 have unblended IRAC [3.6/im] 
and [4.5/im] counterparts. In figures [5]|8] we show the 

^ http:/ /ssc. spitzer. caltech.edu/irac/dh/ 
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Fig. 4. — Limiting stellar mass for galaxies with a stellar pop- 
ulation formed in a burst at 2; = 7 with observed magnitude 
HiQo = 26.0 or Kg = 23.2. The HUDF data are sensitive to pas- 
sively evolving stellar populations with masses above these limits. 

ACS, NICMOS, and IRAC images and spectral energy 
distributions (SEDs) of these objects. 

Including the IRAC data extends our wavelength cov- 
erage, allowing us to discriminate broadly between galax- 
ies dominated by dust-obscured starbursts and galaxies 
with a substantial amount of stellar mass in later-type 
stars (see, e.g., Labbe et al. 2005). In figure [3l we show a 
^160 — [4.5/im] versus ^775 — i^ieo plot to study the stel- 
lar populations and dust extinction properties for the 
16 galaxies with IRAC photometry. We also plot fidu- 
cial colors of a passively evolving stellar population and 
dust-enshrouded starbursts (with E{B — V) = 0.2 and 
0.6, using the extinction law from Calzetti et al. 2000) 
as a function of stellar population age at z = 2.7. These 
simple models bound the range of colors observed in the 
Jiio — HiQo > 1 mag objects, comparable to findings for 
red, J — i^Tg-selected galaxies in other fields (Labbe et al. 
2005; Papovich et al. 2006). Blue Hiqq - [4.5/im] colors 
(< 2 mag) at this redshift (rest-frame 0.4 — 1.2 micron) 
result from substantial older stellar populations, which 
produce a significant Balmer/4000 A break, causing red- 
der ^775 — HiQo colors (and also the red Jno — i^ieo color 
that satisfy the selection criteria). Objects at these red- 
shifts with redder i^ieo — [4.5/im] colors (> 2 mag) require 
dust-enshrouded, ongoing star formation. The galaxies 
with Jiio — HiQo > 1 and IRAC photometry all have 
a similar range of ^775 — i^ieo colors, making it difficult 
to differentiate between the two types of stellar popu- 
lations. However, these galaxies span a large range in 
^160 — [4.5/im] colors, from 1-4 mag. Many of the objects 
appear to have colors consistent with dust-enshrouded, 
young stellar populations, while a few sources have colors 
consistent with a substantial amount of their stellar mass 
in older stars. Based on figure 3, there are few galaxies 
whose colors are consistent with a pure passively evolv- 
ing, older stellar population. We discuss this further in 

§4. 

We also include other multiwavelength data to inter- 
pret the nature of the Jno — i^ieo > 1 niag sample. In 
particular, we use Spitzer /MIPS 24 /im data for this field 
from the catalog of Papovich et al. (2004). These data al- 
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TABLE 1 

Photometric properties of the NICMOS-UDF Jho - Hieo > 1-0 sample 



NIC ID 


R.A. 


Decl. 


Ztot 


Htot 


B - z 


V - z 


I - z 


z - H 


J - H 


[3.6/im] 


[4.5/im] 


Alternate 




[J2000] 


[J2000] 


AB 


AB 


AB 


AB 


AB 


AB 


AB 


AB 


AB 




30* 




-97^d(^^^ 1 S 


23.60 


21.73 


1.55 


1.28 


0.63 


1.87 


1.18 


19.96 


19.77 




66 


no. 00. on 00 


Zi .^U.Z/O.vJ 


25 Al 


24.00 


1.48 


0.64 


0.32 


1.41 


1.10 








86 


no. 09. 07 1 


-Z ( .^u. z 0. y 


24.99 


23.23 


5.41 


2.60 


1.03 


1.76 


1.02 






2^ 


88* 


no. 00. 00 ni 


-97-4fi-9fi S 


25.68 


23.30 


2.27 


1.50 


0.71 


2.37 


1.43 


21.59 


21.44 


3c 


99 


no.oo.o£? 1 7 


-97-4fi-97 7 

Z 1 .^U. Z 1 . 1 


24.58 


23.07 


1.47 


0.80 


0.45 


1.51 


1.05 








111* 


UO.OZ.^-L. ( 


-97-4fi-9S Q 

-Z ( .^U. ZO. c/ 


27.12 


25.53 


1.22 


0.99 


0.24 


1.59 


1.21 


24.44 


24.48 




120* 




-97-4fi-S1 ^ 


26.57 


24.50 


1.53 


0.96 


0.42 


2.07 


1.35 


22.42 


22.21 




215 


no.09.04 77 


-97-4fi-4fi 1 


26.41 


25.08 


1.61 


0.91 


0.23 


1.33 


1.05 








219* 




-97-4fi-4fi Q 


26.16 


23.76 


1.79 


0.76 


0.24 


2.39 


1.78 


22.21 


22.11 


9151^ 


223* 


03 -32 -35 06 


-27-46-47 6 


25.79 


23.36 


1.97 


0.97 


0.43 


2.43 


1.62 


21.61 


21.48 


9024^ 4c 


234 


03:32:42.96 


-27:46:49.9 


23.43 


21.81 


4.29 


2.52 


0.90 


1.62 


1.03 








245 


03:32:39.88 


-27:46:50.0 


26.67 


25.71 


0.07 


0.20 


0.04 


0.97 


1.03 








269* 


03:32:38.73 


-27:46:54.1 


25.74 


24.63 


1.11 


0.45 


0.18 


1.10 


1.09 


23.28 


23.55 




281 


03:32:41.67 


-27:46:55.4 


26.03 


23.72 


2.05 


1.54 


0.97 


2.31 


1.28 








319* 


03:32:34.42 


-27:46:59.7 


24.28 


22.31 


3.03 


2.21 


0.93 


1.97 


1.24 


20.40 


20.25 


5^ 


478 


03:32:42.52 


-27:47:14.3 


26.14 


24.75 


1.77 


0.52 


0.04 


1.39 


1.12 








535 


03:32:34.62 


-27:47:20.9 


25.00 


23.38 


2.38 


1.94 


0.84 


1.61 


1.08 








625* 


03:32:43.46 


-27:47:27.4 


26.63 


24.34 


2.63 


0.82 


0.18 


2.29 


1.77 


23.18 


23.26 


6140^ 


683 


03:32:36.96 


-27:47:27.2 


25.22 


23.23 


2.33 


1.49 


0.62 


1.99 


1.08 








706* 


no . 00. /( 'TO 
Uo:oz:4z.7o 


-z7:47:oo.y 


24.37 


22.74 


1.90 


1.51 


0.75 


1.63 


1.05 


21.28 


21.27 


n 

5256 


902 


03:32:38.10 


-27:47:49.8 


24.47 


22.17 


3.82 


1.80 


0.93 


2.30 


1.10 






3650^ 


921* 


03:32:39.11 


-27:47:51.6 


25.50 


23.11 


3.36 


1.55 


0.86 


2.39 


1.19 


21.92 


21.96 


3574^ 13^ 


935 


03:32:33.26 


-27:47:52.4 


27.23 


24.77 


2.05 


1.97 


0.98 


2.46 


1.33 








989* 


03:32:34.93 


-27:47:56.0 


26.21 


24.76 


1.70 


0.57 


0.35 


1.44 


1.06 


23.37 


23.44 




1078* 


03:32:42.87 


-27:48:09.5 


28.06 


24.54 


2.02 


1.03 


0.62 


3.51 


2.49 


22.03 


21.78 


12182^, 1^ 


1082 


03:32:38.70 


-27:48:10.3 


27.35 


25.84 


1.75 


0.73 


0.09 


1.51 


1.34 








1172* 


03:32:41.74 


-27:48:24.9 


27.75 


24.56 


2.02 


2.80 


0.83 


3.18 


1.39 


22.48 


22.32 




1184 


03:32:38.75 


-27:48:27.1 


26.43 


24.57 


1.84 


0.53 


0.17 


1.86 


1.30 






15^ 


1211* 


03:32:39.16 


-27:48:32.4 


25.78 


23.25 


2.30 


1.19 


0.33 


2.52 


2.02 


22.08 


22.06 


1927*^, 1446^ 


1237* 


03:32:38.73 


-27:48:39.9 


27.07 


24.48 


2.37 


0.97 


0.38 


2.59 


1.95 


22.03 


21.82 


12183^, 2^ 


1245 


03:32:39.58 


-27:48:42.2 


26.21 


25.25 


1.49 


0.59 


0.15 


0.95 


1.07 








1267 


03:32:39.66 


-27:48:50.6 


24.36 


22.38 


2.96 


1.01 


0.22 


1.98 


1.57 






1223^ 



Note. — Objects with unblended IRAC photomet ry are marked with *. References for the alternate ID's listed are (a) iChen M^^^^ (l2QQi) . 
(b) Daddi et al.. (.2005) . (c) T oft e t aL (2005) and (d) lYan et al.l (I2004al) . 



low US to test for the emission from dust associated with 
either star formation or obscured AGN in these galaxies. 
We also include Chandra X-ray data from the catalog 
of Alexander et al. (2003). At the redshifts of interest 
for our sample, the X-ray data are sensitive accretion of 
material onto super massive blackholes, and thus allow us 
to test for the presence of AGN in our sample. 

3. SED FITTING: MASS ESTIMATES FROM A SIMPLE 
PHOTOMETRIC REDSHIFT CODE 

3.1. Fitting method 

To study the properties of our selected objects, we fit 
stellar-population synthesis models to the photometry 
and derive redshifts and stellar masses. We model the 
galaxies as the sum of two stellar components, young and 
old, and consider the effects of reddening. Because the 
star-formation history is constrained generally poorly by 
SED fitting, and the stellar population age and dust ex- 
tinction have strong degeneracies (see, e.g., Papovich et 
al. 2001), we consider only non- negative linear superposi- 
tions of these two models. The combinations of these two 
models encompass the range of plausible star-formation 
histories, and in particular, as we discuss below, they al- 
low us to derive a stringent upper limit on the amount 
of stellar mass in old stellar populations. 

We therefore have four model parameters: the redshift, 
dust extinction, mass in the young stellar component, 
an d the mass in the o l der s tellar component. We use 
the lBruzual fc CharlotI (|2003l ) population synthesis code 



to generate the two stellar population components for 
our SED fitting routine: a 100 Myr constant star for- 
mation (composite stellar population: CSP) model and 
a maximally old passively evolving (simple stellar popu- 
lation: SSP) model formed at z = 7, both generated at 
solar metallicity and with a Salpeter initial mass func- 
tion (IMF). Our maximally old model age is a function 
of redshift: the age at a given redshift is equal to the age 
of the universe at that redshift minus 0.75 Gyr, the age 
of the universe at z = 7 . We attenuate the CSP model 
using the Calze tti et al.l ([200 0) dust law; the old model 
is not attenuated. At each z and E( B — V), we compute 
model photometry by convolving the iBruzual fc CharlotI 
(2003) spectra with the 8 filter response curves (e.g., 
Papovich et al."2001), taking in to account the HI opac- 
ity of the intergala ctic medium (IGM) along the line of 
sigh t (iMadaul 119951) . according the models described in 
iFaj] ^999"). 

We fit the object fluxes to this four-parameter model 
assuming that the errors in the fluxes are symmetric 
and Gaussian. In addition to the photometric uncer- 
tainties on the magnitude of each bandpass, we add a 
systematic uncertainty proportional to the flux density 
of cTsys//^ = 0.04, 0.08, and 0.15 for the ACS, NICMOS, 
and IRAC photometry, respectively. These errors repre- 
sent systematic uncertainties in our matched-band pho- 
tometry, as well as the fact that our model does not sam- 
ple the full range of possible parameters infinitesimally. 
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TABLE 2 

Best- FIT parameters with HST data only 



NIC ID 



'^phot 



30* 

66 

86 

88* 

99 

111* 

120* 

215 

219* 

223* 

234 

245 

269* 

281 

319* 

478 

535 

625* 

683 

706* 

902 

921* 

935 

989* 

1078* 

1082 

1172* 

1184 

1211* 

1237* 

1245 

1267 



1.4^ 



/,+0.2 
-0.1 

2.6±0.2 

1 1+3.3 

-'-•-'--0.1 

9 9+0.1 
^•^-0.2 

2.1±0.1 

7+0-1 

^•'-0.0 

2 4+0.3 

0.6±0.0 

2.8±0.1 

2 9+0-3 
Q 9+0.1 

2.4±0.1 

9 7+0.3 
^- ' -0.1 

1 Q + 0.2 
-•-•^-0.1 

l.liO.l 
3.1±0.2 
l.liO.l 
3.3±0.4 
1 6+0-2 

1 2+0-1 
1.9±0.1 

2 0+0-1 

1.9±0.2 

3.1±0.3 
2 Q+O-3 

9 7+0.0 
-0.1 

2.0±0.4 
3.2^ 
2.9^ 
3.3^ 
2.8^ 
3.2^ 



n> + 0.2 
-0.3 
□ + 0.4 
-0.2 
^ + 0.3 
-0.4 
3 + 0.3 
-0.1 
n> + 0.2 
-0.3 



E(B-V)" 


log Mo/Mq^ 


log My/MQ^ 




v2 , ^ 
■^u,best 




0.2±0.1 


11.2±0.1 


9.0+0-^ 

— O.o 


QQQ + O-004 
— 0.010 


5.2 


0.0 


0.3+0-3, 

— 0.0 


10.6+O-i 

— 0.4 


9.6-^°? 

— 0.1 


0.915-^°°2i 

— O.zUy 


2.5 


0.3 


1 1+0.3 


in + 1-5 

10.3Io 2 


n Q+0.3 

9-8-0.5 


n '7«Q+0.224 


0.0 


0.8 


U.DZlZU. 1 


in s+o-i 


1 n n-i-n a 


n Qfi7+0.083 
U.8b/'_Q 231 


n 7 


n n 
u.u 


n zi-i-n n 


1 n 7-i-n 1 

lU. ( ZCU. 1 


1 n -1-0 1 

lU. ZCU. 1 


n Qp;7+0-018 
U.50^_Q Q23 


O.U 


n n 
u.u 


n n+o-i 


o.oztU.U 


fi p;+0.3 


n qqfi+0-000 
^•^^'-'-0.005 


1 7 q 


n n 
u.u 


n A -i-o 1 

U.^tztU. i 


1^-^-0.1 


Q 9 + 0.3 


n q48+0-032 


n A 


n n 
u.u 


U.OZCU.U 


a fi+0.3 




n n4S+o-05i 


1 9 


n n 
u.u 


n Q_|_n n 


1 1 D-i-n 1 

ll.UZlZU. 1 


q Q+0.0 


n q«i +0-004 
^•^^-^-0.002 


0.0 


n n 
u.u 


n Q_|_n 1 

U.OZlZU. 1 


n -^+0-2 


q 4+0-2 


n q«fi+0-009 
U.y»b_Q Q27 


1 

1.0 


1 1 
1 . 1 


1 1+0-2 


lU.DZlZU. 1 


1 n 1 -l-n 9 
lU. izcu.z 


PI 700+0.094 


1 R 
1.0 


n n 
u.u 


0±0 


10 ±0 


8.1±0.0 


989±0 001 


14.7 


0.0 


0.2+g-; 


10.4-^°-? 

— U. i 


9.2+°-° 


0.941+°-°^3 

— U.U2U 


2.3 


0.0 


0.3+0-? 

— 0.1 


10.7±0.1 


8^7+°-! 

— 0.5 


QQ2 + 0-005 
— 0.029 


0.6 


0.0 


— 0.0 


10.1-^°-? 

— 0.5 


10.2±0.2 


0.421-^°?i« 

— 0.200 


1.9 


0.0 


0.2+°l, 

— u.u 


10.7-^03 


9-2+0? 


— U.UOD 


3.1 


2.5 


— 0.5 


io.o-^°-3 

— 0.5 


9.3+?-3 

— 1.5 


n 7QQ+0.208 

U. 1 00 r> oon 

— 0.330 


4.4 


0.0 


0.2±0.1 




8.9±0.4 


0.993+°™? 


3.9 


0.0 


0.6±0.1 


— U.o 


— U.O 


0.800-^°i39 

— U. zy i 


0.6 


0.0 


^^ — 0.2 


10.6±0.1 


' —0.6 


988+°-°2? 


3.5 


0.7 


0.8±0.1 


— O.o 


10.8-^°-2 

— 0.4 


0.682+noQ? 
— 0.295 


2.8 


0.0 


0.7±0.1 


10.8-^°-^ 

— O.o 


10.2+°i 

— 0.5 


O.oll rt ona 
— 0.3Uo 


1.3 


0.0 


0.5±0.4 


10.3l°-J 


8.6±1.2 


0.9781^45 


0.9 


0.7 


0.2±0.1 


10 7+0-2 

' -0.3 


9.1±0.3 


979+0-014 

U.»^»_Q Q4g 


1.6 


0.1 


O.liO.l 


10.8±0.2 


7 c:+0.6 
' -^-0.4 


1 000+0-000 

i.UUU_Q QQJ^ 


6.6 


0.6 


0.4±0.0 


Q 6 + 0.1 
y.o_Q 2 


9.2±0.0 


707+0-062 
^- -0.112 


6.0 


0.0 


1 1+0.3 

-•-•-•--o.s 


10 4+0-2 


io.il§:^3 


c^-, 0+0. 485 


1.0 


0.8 


0.2±0.1 


10.9±0.2 


9.0±0.3 


n qS7+0-009 
^.i^oi _o.029 


0.4 


0.2 


0.4±0.1 


11 2+0-3 

-•--•-•^-0.2 


9.8±0.3 


960+0-028 
u.ybu_Q Qgj^ 


9.5 


1.0 


1+0-2 


11.0±0.2 


O q + 0.5 


n qqs+O-OOi 
0.yy8_o.oo8 


2.7 


0.0 


Q o + O.O 


9.8±0.5 




766+0-182 
U. /bb_Q 260 


5.6 


1.7 


0.4±0.1 


11 6+0-2 


10.6±0.3 


921+0-054 

^•^^-^-0.180 


5.1 


0.0 



Note. — Columns 2 through 4 list the median values and the 16% and 82% intervals for the fitted parameters. Objects with unblended IRAC 
photometry are marked with an asterisk symbol. Objects 30, 88, 1237 and 1267 are detected in x-rays. 
^ The E(B-V) extinction is applied only to the young component 
^ The model fit mass in the old (SSP) component. 
*^ The model fit mass in the young (constant SFR) component. 
^ The fraction by mass of the old component relative to the total fitted mass. 
^ The best-fit model xt listed. 

^ The difference between the value for the model listed in columns 2 through 4 and the best x^'- 
Ax^ = xliedian- Xbest (also See Section 3.1). 



These systematic errors are added in quadrature to ei- 
ther the fitted errors, for the ACS photometry, or the 
catalog errors for the NICMOS and IRAC photometry. 
We use the resulting statistic to produce the likeli- 
hood C (X exp(— and we take the likelihood to be 
the probability distribution for the 4 parameters given 
the observed fluxes. 

We derive mean values and variances for the 4 pa- 
rameters (and functions thereof) by integrating over the 
full 4-dimensional likelihood. The integrations are done 
with a hybrid grid plus Monte Carlo approach. The 
is quadratic in two parameters (the old and young 
masses, Mq and My) and more complicated in the other 
two (redshift and reddening). We construct a grid in 
redshift and reddening over the range < z < 5 and 
< E{B — V) < 1.5. For each point, we use linear 
methods to find the best fit and Gaussian covariance ma- 



trix for the masses. We generate Gaussian-distributed 
points given this covariance region and weight points by 
exp(— Xbest /^) best-fit x^, dropping any points 

that have negative masses. The concatenation of all of 
the points is a Monte Carlo of the full 4-dimensional 
likelihood, and any function can be computed by a sim- 
ple weighted sum of the function evaluated at all of the 
points. Similarly, we can construct quantiles by sort- 
ing the points according to a given function evaluation 
and cumulating the weights. By using a grid in the non- 
quadratic directions, we ensure that multiple minima are 
found and correctly weighted. 

The derived physical parameters are presented in ta- 
bleland table [3l In these tables we define /om to be the 
ratio of the mass in the passively evolving component. 
Mo to the total mass, Mq + My, where My is the mass 
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TABLE 3 

Best-fit parameters for objects with IRAC photometry 



NIC ID 


'^phot 


E(B-V)" 


log Mo/Mo^ 


log M^/Mo^ 




v2 , ^ 




30** 


1-8+0? 


0.3+g-J 


11.5±0.1 


9.7+0-^ 


0-984+^26 


9.2 


0.0 


88** 


2.2+°? 

— 0.1 


0.7+0-? 

— 0.1 


10.8±0.2 


10.3+°-^ 

— 0.2 


0.718+O-J,^? 

— 0.071 


0.7 


0.0 


111 


— u.u 


0.4±0.0 


9.6±0.1 


8.5+°-i 

— 0.0 


0.914+°°;° 

— 0.015 


11.3 


0.0 


120** 




0.2±0.1 


ll.liO.l 


8.9±0.4 


0.994+°°?^ 

— 0.012 


2.0 


0.2 


219* 




0.3+n 


ll.liO.l 


Q Q + o.o 

3 


982+°-°ii 

U.o/O.^ Q QQ2 


2.2 


0.0 


223* 


S.litO.l 


0.2+0-1 

— 0.0 


11.4±0.1 


9.1+n'n 

— 0.0 


0.995+0-000 

— 0.007 


0.7 


0.0 


269* 




0.2±0.0 


10.4±0.1 


9.2±0.0 


944+0-007 

^•^^^-0.008 


1.9 


0.0 


319** 


1 9+0.0 


^•^-0.0 




10.3dz0.1 


Q 044+0.129 


1.3 


0.0 


625* 


2.7±0.1 


4+0-0 


10.6±0.1 


9 4+°-i 


Q07+0.039 

-0.015 


3.3 


0.0 


706 


1 q + 0.2 


0.6±0.0 


10 3+°-2 


9.9±0.1 


703+0-095 

^- '^"^-0.137 


4.2 


0.0 


921 


1.9±0.1 


0.5±0.1 


10.9±0.1 


9.5±0.3 


966+°-°^^ 


1.6 


0.0 


989* 


2 9+°-^ 


Q o + O.O 


10.5±°:? 


Q 3 + 0.1 


qop: + 0.042 


1.3 


1.5 


1078 


2 9+°-'^ 

^•^-0.3 




io.6±°:^3 




0.386±0-?i4 


5.8 


2.4 


1172 


9 1+0.2 


I 1 + 0.3 


10 4+°-i 


io.i±°-i 


661+°-°'^^ 

U.ODl_o.o72 


1.0 


0.3 


1211* 


2.7±0.1 




ll.OiO.l 


9.9±0.2 


949+0-013 

"^•^^^-0.073 


5.0 


0.5 


1237** 


o ^ + 0.2 
"^•"-0.0 


o+°-i 


11 2+°-i 


8 0+°-^ 


QQQ + 0-000 
^•^^^-0.001 


1.9 


0.0 



Note. — Columns 2 through 4 list the median values and the 16% and 82% intervals for the fitted parameters. The six galaxies with broad- 
band SEDs consistent with distant evolved stellar populations are marked with a star symbol. Objects detected at 24/im are indicated with a 
double-asterisk symbol. Objects 30, 88, and 1237 are detected in x-rays. 
^ The E(B-V) extinction is applied only to the young component 
^ The model fit mass in the old (SSP) component. 
*^ The model fit mass in the young (constant SFR) component. 
^ The fraction by mass of the old component relative to the total fitted mass. 
^ The best-fit model xt is listed. 

^ The difference between the value for the model listed in columns 2 through 4 and the best x^- 
Ax^ = xLed^ar^- xLst (also See Sectiou 3.1). 



in the star- forming component, 

/oid = Mo/(Mo + My). (1) 
We characterize the probabihty distribution function for 
the various parameters by the median and the 16% and 
84% intervals (equivalent to the la intervals for a Gaus- 
sian distribution). However, in cases with multiple min- 
ima these simple quantiles will not fully represent the 
distribution. 

3.2. Fitting results 

We apply this SED fitting method to derive physical 
parameters for our objects. We present the results of the 
ACS and NICMOS photometry fitting analysis in detail 
in Tableland summarize them here as follows. 

The fits to 11 objects in our sample favor models with 
relatively lower redshifts (^hot < 2.5) and lower stellar 
mass fractions in old stellar populations (/om < 0.9): 
NIC IDs 86, 88, 99, 215, 234, 319, 535, 683, 902, 921, and 
1172. We note that on average this set of galaxies has 
more dust, with a mean fitted E{B — V) ~ 0.8, than the 
entire sample, which has a mean fitted E{B — V) ^ 0.4. 
The fits to 7 objects favor models with lower redshifts 
and higher mass fractions in older stars (^^phot < 2.5 and 
/old > 0.9): 30, 111, 120, 245, 281, 706, and 935. For 
objects 1082 and 1245, the models favor a solution with 
high redshift (%hot > 2.5) and a relatively low fraction 
of mass in old stars (/om < 0.9). For 12 objects the 
models favor solutions with high redshift (%hot > 2.5) 
and a relatively high fraction of stars in the old stellar 
population (/om > 0.9): 66, 219, 223, 269, 478, 625, 989, 
1078, 1184, 1211, 1237 and 1267. 

We find that when we include the IRAC [3.6/im] and 
[4.5/im] data in the model fitting, 12 objects (of the 16 



with unblended IRAC photometry) have consistent re- 
sults compared to the model fits without the IRAC data. 
Of these 12 objects, 6 objects are fitted with high-redshift 
model galaxies with only small amounts of rest-frame 
UV light associated with a young stellar population or 
on-going star formation (for this set 2^phot 

> 2.7 and 

/old > 0.93). The fitted model parameters for these 16 
objects are listed in Table [3l Here we discuss the particu- 
larities of these objects and their corresponding fits. We 
will not individually discuss the objects without IRAC 
photometry. However, we note that object 1267, is de- 
tected in the X-ray data, with a 0.5 — 2.0 keV flux of 
fx = 7.5 X 10"-^^ ergs cm~^ s~-^. 

NIC ID = 30: The difference in the fitted param- 
eters when we include the IRAC data is significant: 
^Zphot ^ 0.4. As expected, the masses change as well. 
However, both values for /qm remain high. The best-fit 
is very large for the IRAC fit (x^ = 9.2). This is due 
to the large amount of IRAC flux that our models inade- 
quately characterize (see figure[5j); this excess i s also driv- 
ing the selection of a higher redshift and mass. iToft et al] 
(2005) derive a photometric redshift of z = 2.2, roughly 
consistent with our IRAC fit redshift z = 1.8lo"^. How- 
ever, this object is detected with MIPS at 24/im with 
/^y (24/im) = 190 /iJy, and it has an X-ray data detec- 
tion of fx {0.5 - 2keV) = 8.4 x 10"^^ ergs cm-^ s-^ 
Based on this data it is likely that this object likely has 
an AGN; therefore we consider the lower redshift fit more 
plausible. Furthermore, it is possible that the AGN con- 
tributes to the broadband photometry, and our model 
fits may not be valid for this object. 
NIC ID = 88: The two fits, with and without the IRAC 
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NIC ID=30 



^110 



[3.6^tm] 



[4.5^tm] 



20 r 



m 22 
< 

^ 24 
I 26 
28 
30 
32 



I ' — ' — ' — ' — ' — r~ 

r • Doto 
: xBest fit 

i- 

Hi — ^-'•■^ 


r- 


^ 


— 1 1 1 1 — 1 — I 












^phol ~ 




- 1 




froc = 


n QBA +0.002 - 
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Fig. 5. — Gallery of Jno — Hiqq > 1.0 and i^ieo < 26.0 selected objects; we show imaging, broad-band SED's and the best fit model 
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IRAC imaging is also displayed at constant fu per arcsec on a log scale but on a different relative scale than the S435 V606^775^850-^ii0-f^i60 
imaging. The broad-band SED's are shown in the lower panels with solid circles, la error bars are indicated. The crosses indicate the best 
fit model. The dashed line indicates the fraction, by mass, of light that is attributed by the best fit to the young model; similarly, the solid 
line indicates the fraction attributed to the old model. 



data, are consistent, with Zphot ~ 2.2 and /om ~ 0.8. It 
is detected at 24/im with /j,(24/im) = 70 /iJy, although 
inspection of the 24/im image shows that the object is 
blended (with a PSF FWHM 6''). This object also 
has a 0.5 - 2.0 keV x-ray flux f x = 8.8 x 10"^^ ergs 
cm~^ s~^ (jAlexander et al.lf2003[ ). In conjunction these 
data indicate that this object is likely host to an AGN, a 
conclusion that is supported by the fact that this object 
has red IRAC [5.8/im] — [8.0/im] colors and that it has a 
compact morphology in the HST images (see flgure[5|). 
As with NIC ID 30, the presence of an AGN may imply 
that the models are inadequate. 

NIC ID = 111: The best-flt model, derived including 
the IRAC data, is lower redshift older stellar population: 
%hot = 1.5 and /om = 0.914. The model derived from 



the ACS and NICMOS data alone have a similar /qm, but 
the redshift is much lower, z = 0.7. This discrepancy is 
driven by the bright IRAC photometry; our models are 
simply not red enough at z = 0.7 to reproduce the mid-IR 
data. This indicates that the ACS and NICMOS broad- 
band photometry alone does not adequately characterize 
this object. 

NIC ID = 120: The fltted photometric redshifts and 
masses are consistent with a distant old stellar popula- 
tion dominating the flux in the rest-frame optical and a 
younger stellar population (or on-going star formation) 
dominating in the rest-frame UV, with %hot = 3.4 and 
/old = 0.994. The physical parameters derived with- 
out the IRAC data are not consistent with this inter- 
pretation, with a much smaller distance {z = 2.4) and 
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a somewhat lower fraction (/qm = 0.948). Inspection of 
the images reveals that the blue component dominating 
the rest-frame UV photometry is unresolved and just off- 
center (see figure [5j). Therefore, the possibility that this 
SED is the result of the projection of more than one ob- 
ject cannot be ruled out. Despite the fact that [5.8/im] 
and [8.0/im] IRAC fiuxes are in good agreement with the 
distant old population scenario, we do not include this 
object in our set of distant and evolved object candidates 
because of the discrepancy between the two model fits. 
Furthermore, this object is detected at 24/im. 
NIC ID = 219: Our fitted physical parameters are con- 
sistent with a high redshift galaxy dominated by an old 
stellar population (%hot = 2.8 and /om = 0.98). The fit 
derived from the dataset including IRAC data is consis- 
tent with this model. The derived redshift is also con- 
sistent with the photometric redshifts of Chen & Marzke 
(|2004f ). 

NIC ID = 223: The fitted model is a distant galaxy 



dominated by an old stellar population, with best-fit pa- 
rameters 2;phot = 3.1 ± 0.1 and /om = 0.995, consistent 
with t he model derived from t he ACS and NICMOS data 
alone. I Chen fc Marzkd (|2004 l) derive a photometric red- 
shift of 2;p hot = 3.11,^ i n ^ood agreement with our derived 
redshift. Yan et al. (2 004al ) derive a 2;phot = 2.9 and 
Toft et^l. (2005) derive a 2;phot = 3.4, however, these 
authors do not quote errors, making it hard to compare 
with our model. 

NIC ID = 269: The fit parameters are consistent with 
a distant object dominated by an old stellar population 
with a small amount of ongoing star formation, for both 
fits with and without IRAC data (%hot = 2.7 and /om = 
0.94) 

NIC ID = 319: The most-likely model is a low red- 
shift dusty starb urst: 2;phot = 1-2, E{B — V) = 0.9 
and /old = 0.34. iToft et all (120051 ) derive a Zphot = 1-8 
for this object. This object is detected at 24/im with 
/^y (24/im) = 110 /iJy, supporting the interpretation that 



10 




INIC 10=706 




[3.6^tm] 



[4.5Atm] 



I ' ' ' ' ' — 1 ' 

r • Doto 
: xBest fit 


1 1 1 1 1 1 1 — ; 








z = 1 3-'°-2 -. 
froc = 0.703:g;095 "i 



2 

A (/xm) 



NIC 10=989 



[3.6^tm] 



[4.5/zm] 




















[3.6^trT 






[4.5/im] 



20 r •Doto 

xBest fit 



— 1 1 1— 



— 1 1 1 1— 




+0.1 

froc = 0.935!g;0t2 



2 

X (/xm) 



Fig. 7. — Same as figure [5] 



this object is a lower redshift dust-enshrouded starburst. 
NIC ID = 625: The most-hkely model is consistent with 
a distant galaxy dominated by mass by an old stellar pop- 
ulation with a small amount of ongoing star formation 
{zphot = 2.7 ± 0.1 and /om = 0.94). The redshift derived 
without the IRAC data is higher, although poorly con- 
strained, 2;phot = 3.3 ±0.4, and the stell ar population re- 
mains dominated by an old population. IChen fc Marzkd 
(|2004l ) derive a redshift for this object of ^hot = 3.33, 
consistent with our fit excluding the IRAC data (see ta- 
ble 2). 

NIC ID — 706: The most-likely model i s consistent with 
a low- redshift dusty star forming galaxy. I Chen fc Marzkd 
(|2004l ) reach a similar conclusion for this object. 
NIC ID = 921: The fitting results are consistent with an 
intermediate redshift old stellar population, with %hot = 
1.9 ± 0.1 and /om = 0.97. Daddi et al. (2005) measure 
a co nsistent result with a spectroscopic redshift of z = 
1.98. lYan et aP (|2004af ) also derive a similar photometric 



redshift. Our mass, 8 x 10-^^ M^:^ is i n good agreement 
with the mass derived bv iDaddi et al.l ((2005). 
NIC ID = 989: The resulting model fit is consistent 
with a distant galaxy dominated by mass by an old stel- 
lar population and with a small amount of ongoing star 
formation (%hot = 2.9 and /om = 0.94). Exclusion of the 
IRAC data yields a similar model. 

NIC ID = 1078: The most-likely model parameters de- 
rived from the HST+IRAC data are consistent with a 
distant, 2;phot = 2.91q"3, Y^^S dusty stellar popu- 
lation, with a fitted /om = 0.4 and E{B - V) = 1.0, al- 
though the redshift is not particularly well constrained. 
Fitting the photometry while excluding the IRAC data 
gives a consistent redshift, but a very different interpreta- 
tion of this object. In this case the model is dominated by 
an old stellar population, with /qm = 1-0. The relatively 
bright IRAC flux densities, taken in conjunction with 
the HST photometry, are inconsistent with our models of 
old stellar populations, which require a flattening of the 
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SEP at wavelengths blueward of about 0.5/im rest-frame. 
IChen fc M arzke (2004) derive a 2;phot = 2.91, in good 
agreement with our results. However, they find that this 
object is consistent with a E/SO SED, far more evolved 
than the constant star- f orming model we d erive when we 
include the IRAC data. lYan et all ('2004a') select this ob- 
ject using a fu{S.6jiim)/fjy{zs5o) > 20 selection and derive 
a %hot = 3.6, also cons istent with our poorly constrained 
Zphot- lYan et aP (|2004a) do not include this object in 
their analysis due to the large uncert ainties in their fit- 
ted redshift. However, we note that Yan et all (|2004al ) 
report a non-detection for this object (2^850 > 29.5 mag; 
using data from Yan & Windhorst 2004c), while we re- 
port a detection of ^850= 28.1 mag. Upon inspection 
of the unconvolved ACS z-band image, we see low sur- 
face brightness fiux at the location of this object and 
we measure a detectable fiux in our PSF-convolved im- 
age. Incidentally, although the analysis of Mobasher et 
al. (2005) supports a very high redshift interpretation for 



this object (z > 6), subsequent analyses by Dunlop et al. 
(2007, MNRAS, 376, 1054) and Chary et al. (2007) favor 
an interpretation that this source resides at substantially 
lower redshift z < 4, consistent with our interpretation. 
NIC ID = 1172: This objects' most-likely model pa- 
rameters are consistent with an intermediate redshift 
old stellar population with a significant amount of UV- 
bright younger st ars: ^hot = 2.1^3^ and /om = 0.66. 



2.7, somewhat larger 



lYan et all (|2Q04aD derive a Zph, 
than our value of 2;phot- 
NIC ID = 1211: The resulting model fit is consistent 
with a distant galaxy whose mass is dominated by an 
old stellar population. We derive a ^^nho t = 2.7 =b 0.1 
and an /qih = 0.95 . [Chen Marzkd (l2004l ). |paddi et all 
(|2005r ). iToft et all (|2005r ) and lYan et all |2004a) derive 
Zphot = 3.43,2.47,3.24,2.8 respectively. We suspect the 
discrepancies between the different fits result from as- 
sumptions about the star-formation histories of the mod- 
els and differences in the photometry. Our derived mass. 
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Mtot = 1-1 X 10^^ Mq, is highe r than, but not inc onsis- 
tent with, the mass derived by iDaddi et alJ (l2QQ5f). Al - 
though this object is selected bv lChen fc Marzkd (|2QQ4l ) 
they do not report its mass as they classify it as a 
dusty starburst. This object has an upper limit X-ray 
flux of / (0.5 - 2keV) = 3.7 x 10"^^ ergs cm-^ s"^ 
([Alexander et al.ll2QQ3l ), implying the possible presence 
of an AGN. However, the rest-frame UV through IR pho- 
tometry is consistent with with integrated stellar popu- 
lations and therefore we include this object in our sample 
of distant, evolved galaxies. 

NIC ID = 1237: This objects' best-flt model is consis- 
tent with a distant, old stellar population: Zphot = 3.6 
and /old = 0.999. The exclusion if the I RAQ data gives a 
similar flt, ^phot = 3.3 and /om = 0.998. IChen Mar^ 
( [2004 ) derive a much higher red shift for this object, 
Zphot = 4.25. lYan et alT (|2004at ) also select this ob- 
ject but do not consider it due to the fact that their 
%hot = 3.4 is very uncertain. We note that our photom- 
etry for this object is likely biased because of the prox- 
imity of a very bright spiral galaxy, just visible in the 
lower left corner of flg. 8. It is likely that the image con- 
volution causes low-surface-brightness emission for this 
objects to enter our aperture. Furthermore, our pho- 
tometry could be biased by the presence of nearby, faint 
sources in the ACS data that might be blended with this 
source at the NICMOS H-band resolution. This object 
also has a marginal 24/im detection due to the fact that 
it is blended, with /^^ (24/im) = 76 /iJy, which may call 
in to question the validity of the best-flt model. Because 
of these reasons we reject this object as a candidate for 
a distant, evolved galaxy. 

We note that although our photometry is generally in 
good agreement with other published values which we 
are aware of, a few of our objects have signiflcant differ- 
ences (notably NIC ID 1078 and 1237). In the case of 
NIC ID 1078, we suspect our photometry is appropriate 
for our NICMOS-selected catalog. For NIC 1237, scat- 
tered light from a nearby, large galaxy likely affects our 
photometry measurements. Regardless, neither of these 
objects enters our flnal sample for galaxies whose mass is 
dominated by older stellar populations, so these do not 
affect our flnal conclusions. 

We derive masses assuming a Salpeter IMF, where the 
IBruzual fc CharlotI (|2003[ ) models use upper and lower 
mass cutoffs of lOOJVlq and 0.1 M©. For an SSP 
IBruzual fc CharlotI (|2003[ ) model, assuming a Chabrier 
IMF would decrease our masses by a factor of 1.8, rel- 
ative a Salpeter IMF and for the same assumed upper 
and lower mass cutoffs listed above. For a more detailed 
discussion of how masses depend on the IM F parameters 
such as the slope and the mass cutoffs, see P apovich et al.l 
(2001). We note that Chen fc Marzke (200j) identify red 
objects in the HUDF and derive masses for six early- type 
candidates, 4 of which are systematically lower than the 
masses derived here. Our overlap objects are 219, 223, 
625, 1078, and 1237. The total H-band magnitudes also 
differ for these objects; as expected, there is a rough 
trend between increasing flux discrepancy and increasing 
mass discrepancies. Object 223 has the largest difference 
in mass, a factor of ~ 40. Objects 219, 625 and 1078 have 
masses that di ffer by f actors of 18.0, 2.1 and 1.4, respec- 
tively, from thel CheiTfc Marzke (2004) masses, although 
the two estimates for objects 625 and 1078 are entirely 



consistent. Obj ect 1237 is the only for which our mass is 
lower than the iChen fc Marzkd (|2004l ) mass, by a factor 
of ~ 1.2, but as stated previously, this likely results from 
the higher redshift derived by Chen fc Marzke. 

3.3. Other models 

To investigate possible systematic dependences due to 
our modeling, we explore the effects of changing our 
model assumptions. 

We flrst consider a single component flt using the same 
young model described above, a constant star forming 
model attenuated with the iCalzetti et al.i (^2000 ) law. 
This model does not characterize our objects well, yield- 
ing very large values compared to the two-component 
base model, with an average increase in ~ 55. 

We also test a m ore complex dust model, the 
IWitt fc GordonI (|2000f ) SMC shell geometry dust model. 
We try both a single component and a two-component 
flt with this dust model. The single component flt, 
a 100 Myr old constant star forming model with dust 
as a free parameter, does not characterize the data 
well; the best-flt values are very large compared to 
our base mod e l and to the single component flt using 
ICalzetti et al.l (|2000f ) dust discussed aloove. The two 
component fit, the second component being a maximally 
old SSP model, results in an average x^ difference of 

(Xbest, base model " Xbest) = '^'^ O^Cr ah 16 objcctS, 

somewhat worse than our base model. This dust model 
gives lower Xbest ^^l^^s for 9 objects. Object 921 has 
the largest Ax^ of 3.3. This dust model assigns sys- 
tematically higher E{B — V) values, in particular to ob- 
jects 30, 625, 921 and 1211, with best-flt E{B - V) = 
1.4, 2.5, 5.3 and 2.8 respectively. For completeness we 
note that the dust grid over w hich our ni odels span is 
the grid provided by the Wit t fc GordonI (2000) look- 
up tables for this dust law. It is not the same grid as 
that in our base model and it is not a regularly vary- 
ing grid, having values of E(yB — V) ranging from 0.0 
to 12.3, in steps of /\E{B - V) = 0.0875 at low values 
and in steps of AE{B — V) = 1.75 at high values. This 
difference may hamper our ability to make a fair com- 
parison with our base model. However, we do not think 
this is the case as none of our original model flts have 
large values of E{B — V), all objects having model flt 
values of E{B - V) < 1.1. We note that the best-flt 
masses differ by less than a factor of 3.6 between this 
model and our base model, excluding object 319, whose 
mass increases by about an order of magnitude relative 
to our base model. Similarly, the fractions change by 
small amounts: (|A/oid|) = 0.08 over the set of 16 ob- 
jects. The only objects with a signiflcant change in red- 
shift are 88, 319 and 1172 with a Azphot = 0.9, 1.1 and 
1.1 respectively. None of these objects are included in 
our flnal sample of galaxies whose mass is dominated by 
older stellar populations. 

We next consider the possibility that our objects have 
metallicities lower than solar. Again, we perform a sim- 
ilar analysis as outlined above, with both the young and 
old mode ls at the same a^es a s for our base model and 
using the ICalzetti et al.l (|2000l ) dust law, but with both 
components at a metallicity of 0.2Zq. For this model, 
the average difference in x^ is (xLst, base model "Xbest) = 
—2.1, i.e., our base model does a better job at flt- 



13 



TABLE 4 

Comparison between Table 3 and the t-model fitting 



NIC ID 


'^phot 


phot 


log Mtot/Mq 


log Mtot/Mq-^ 


^old 


Age (Gyr)^ 


r (Gyr) 


"^v>,hest 


2 r 


30 


1.8 


1.5±0.1 


11.5 


11 2±0 1 


0.984 


1 47±0 36 


1 19+°-Ii 


9.2 


8.4 


88 


2.2 




10.9 


10.8l°i 


0.718 


omtrot 


^•^^-0.01 


0.7 


1.3 


111 


1.5 


0.8±0.1 


9.6 


7.81°:° 


0.914 


0.14±0.11 


0.13±0.10 


11.3 


4.5 


120 


3.4 


2.3±0.2 


11.1 


10.5l°-j 


0.994 


QQ+0.92 


Ql+9.92 
^•^-^-0.90 


2.0 


1.2 


219* 


2.8 


2.5±0.1 


11.1 


10.9±0.1 


0.982 




0.36±0.16 


2.2 


2.0 


OOQ* 

zzo 


1 


9 O+0.7 


114 


111 A-Ci 1 
1 1. imu. 1 




p 0^+0.44 


n SI +3-09 


n 7 


1 n 

l.U 


269* 

319 


2.7 
1.2 


2.6±0.1 

1 0+°-° 


10.4 
10.5 


lO.liO.l 
10.4lg-i 


0.944 
0.344 


0.75+0-^^ 

— 0.22 

o-o8ig:g^e 


1.68+i^o 

— 1.13 
U.Z0_Q 25 


1.9 
1.3 


2.2 
0.9 


625* 


2.7 


2.7±0.1 


10.6 


10.3±0.1 


0.937 


lQ+0-23 


04+°-°^ 

U.U^_Q Q3 


3.3 


1.1 


706 


1.3 


1.2±0.1 


10.4 


10.2±0.3 


0.703 


n -,0+0.60 
^•-'-'-•-0.09 


7l+'^-i'^ 


4.2 


2.8 


921 
989* 
1078 
1172 


1.9 
2.9 
2.9 
2.1 


1 Q+O-2 

-L-y_i.4 

2.5±0.2 

9 7+0.2 
-0.1 

1 4+0-3 


10.9 
10.5 
11.0 
10.6 


10.6+?-^ 

10 1+°-^ 

11.2±0.1 
10.4ig-j 


0.966 
0.935 
0.386 
0.661 


n 9(^+0.22 

Q 47+0.34 
U.^' _0.37 

U-O^_0.41 


0.68l^:^e? 
^•-•-^-0.11 


1.6 
1.3 
5.8 
1.0 


1.6 
2.2 
1.2 
0.5 


1211* 


2.7 


2.8±0.1 


11.1 


10.9±0.1 


0.949 


0.34±0.11 


0.09±0.01 


5.0 


1.7 


1237 


3.6 


9 7+0.2 
-0.1 


11.2 


11 2+°-^ 

-'--'-•^-0.2 


0.999 


1 Qo + 0.66 
-'-•^'^-0.41 


72+0-^^ 

^- '^-0.45 


1.9 


0.5 



Note. — The columns with r-model parameters are indicated with a r symbol. The six galaxies with broad-band SEDs consistent with distant 
evolved stellar populations are marked with a star symbol. 



ting our object SEDs. This set of low-metallicity mod- 
els yields a lower Xbest ^^^^ objects. For 
the entire set of 16 objects the change in redshift is 
(|Azphot|) = 0-3, with object 30 having by far the largest 
difference: Zphot, base model - %hot = -1.3. The best- 
fit masses change little, with a variation smaller than a 
factor of 3 for all objects. The stellar populations re- 
main consistent for all 16 objects except 1237, which is 
fitted as a young starburst with /qm = 0.07, as opposed 
to our base model fit with /qm, base model = 0.999. This 
result lends credence to the interpretation supported by 
the 24/im and X-ray data that this object may host an 
obscured AGN. We note that altering the metallicity of 
the base model does not significantly change the fitting 
results for the 6 objects we select as candidates for dis- 
tant, evolved galaxy candidates and therefore will not 
affect our final conclusions. 

We also consider the effects of a dusty old population; 
we fit our data using th e two-compone nt base model out- 
lined above and we use ICalzetti et al.l (2000) dust to at- 
tenuate the old component with E{B — V)= 0.05 (or an 
Ay = 0.2). For this model we calculate an average 

change of (xLt, base modei-xLt) = ^-l- We derive xLt 
that are lower than those derived with our base mod- 
els for 11 objects (NIC ID = 30, 88, 111, 120, 269, 625, 
706, 986, 1078, 1172, and 1211). For this sub-set, objects 
111, 625 1078, and 1211 have the largest differences, with 
Ax^ = 32.9, 9.0, 11.7, and 6.2 respectively. The other 7 
objects all have Ax^ < 3.3. The photometric redshifts 
change little, with a (|Azphot|) = 0.15 for all objects, ex- 
cept 111, which has a Azphot = 1.9. The mass fractions 
in old and young stars change little, with (|A/oid|) = 0.03 
for all objects except 219 and 1237, both of which have 
a decrease in /om of /om, base model - /old = 0.9. For the 
6 objects we consider distant, evolved galaxy candidates, 
the derived masses change by less than 50%. 

Finally, we consider single-component r-models with 
varying ages similar to other studies in the literature. 
We summarize the fitting results here and compare, in 



more detail, the best-fit r-model parameters with our 
base model results in Table 4. We generate models over 
a grid in ages and r values, were r is the exponential 
timescale for decay in star- formation. Our age grid values 
are 0.01, 0.02, 0.04, 0.06, 0.08, 0.10, 0.20, 0.50, 0.70, 1.00, 
and 2.00 Gyr, and our r grid spans 0.001 Gyr to 20.0 Gyr 
in quasi- logarithmic steps. We use the same redshift and 
dust grid as that of our base model and we assuni e the 
same dust attenuation law, the ICalzetti et all (|2000l ) law. 
We fit the observed SEDs of our selected objects using 
this 4-dimensional grid of models. With these models we 
obtain lower Xbest values, relative to our base model, for 
11 objects (NIC ID = 30, 111, 120, 219, 319, 625, 706, 
1078, 1172, 1211, and 1237), with objects 111, 1078, and 
1211 having the largest differences of Ax^ = 27, 18, and 
13, respectively. For all 16 objects we measure a change 

in X^ of (Xbest, base model " Xbest) = ^.0; if WC CXcludc 

objects 111, 1078, and 1211, the change in x^ becomes 

(Xbest, base model " Xbest) = The maSSCS of the 11 

objects with lower Xbest r change by less than a factor 2 
with respect to the base model fit, the only exceptions 
being objects 111 and 120, which have factors of ~60 and 
~3 decrease in mass respectively. The fitted redshifts do 
not change significantly, (2:phot,base model - %hot) = 0.3, 
with objects 111, 120, and 1172 having the largest change 
of A2;phot ^ 0-7, 1.0, and 0.7 respectively. For the entire 
sample of 16 it is interesting to note that 14/16 objects 
have base model masses that are higher than those de- 
rived with the r-models. We consider this to be a con- 
sequence of the fact that our base model is better suited 
to measure an upper limit on the masses of our objects; 
this is in part because the maximally old component of 
this model has a limiting mass to light ratio. Further- 
more, we note that for the 6 objects with distant evolved 
fits, we derive masses which are between a factor of 1.3 
and 2.5 greater than those derived from the r-model fits. 
Therefore we conclude that our adopted two-component 
base model is better suited to place upper limits on the 
old stellar population census in our objects than the r- 
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model. We note that the two objects with higher r-model 
masses are 1078 and 1237, which are greater by a factor 
of 1.6 and 1.1, respectively, still consistent within the 
errors with our base model. 

We conclude that a two-component fit is better able to 
characterize the data than a one-component fit and that, 
although some objects are better fit by other models, in 
general our base model characterizes our objects well. 
We show that the changes in our fitted redshifts, total 
masses and fractions (/old) are not strongly dependent 
on the assumed metallicity and dust properties of our 
models. 

4. DISCUSSION AND SUMMARY 

We have investigated the physical properties of 32 
galaxies selected with (Jno — i^ieo) > 1.0 colors in the 
HUDF using deep, broadband imaging from HST and 
Spitzer. Exactly half of these objects (16/32) have un- 
blended IRAC [3.6//m] and [4.5//m] photometry. Here, 
we focus on the subset of galaxies with IRAC photometry 
because the data span a longer wavelength baseline, al- 
lowing us to separate more cleanly the objects into those 
whose colors are consistent with dust-enshrouded star- 
bursts and those dominated by old stars. Thus, where we 
derive number densities and mass densities, we include 
a factor of 2 correction for incompleteness due to blend- 
ing in the IRAC image. Furthermore, including both 
the Spitzer /IRAC data allows us to better constrain the 
properties of the galaxies' stellar populations (see fur- 
ther discussion in, e.g., Labbe et al. 2005; Papovich et 
al. 2006). 

4.1. The Properties of Galaxies with Evolved Stellar 
Populations at z ^2.5 

Based on the full SED fitting, we split the subsample 
of objects with IRAC photometry into two broad cat- 
egories: those objects whose broad-band SEDs are well 
represented by our two component model with a substan- 
tial component of evolved stars combined with a young 
star- forming component, and those objects whose inter- 
pretation changes substantially by including the IRAC 
photometry compared to the model fits to only the ACS 
and NICMOS data. For 12/16 (75%) of the objects, in- 
cluding the IRAC data give roughly consistent best-fit 
models as to the ACS and NICMOS data only (NIC IDs 
88, 219, 223, 269, 319, 625, 706, 921, 989, 1172, 1211, 
and 1237). Because we are primarily interested in placing 
constraints on the number density of galaxies dominated 
by passively evolving stellar populations, we focus here 
on the subset of galaxies whose model fits have 2: > 2.5 
and /old > 0.9. That is, the models favor a scenario 
where the vast majority of the galaxy's stellar mass re- 
sides in an older, passively evolving component. There 
are 6 objects satisfying these conditions (NIC IDs 219, 
223, 269, 625, 989, and 1211). 

For three of these objects (NIC IDs 223, 625, 989), we 
have higher-confidence redshifts due to the presence of a 
Lyman "break" between in the ^435 — Vqoq colors. This 
occurs for galaxies with relatively blue rest-frame UV col- 
ors, and a characteristic "break" in the colors as neutral 
hydrogen attenuation from the IGM absorbs the intrin- 
sic fiux shortward of 1216 A (e.g., Madau 1995). For the 
ACS filters, absorption owing to IGM neutral hydrogen 
affects the 5435-band at 2: > 3, and the presence of this 



break drives the estimate of the redshifts for these galax- 
ies (2;phot — 3.5, 3.2, and 3.1, respectively). Interestingly, 
at 2: 3 — 3.5 the Balmer/4000 A is moving through the 
i^ieo-band, diminishing the J — H color. Therefore the 
fact that these galaxies have Jiio — ^i60 ^ 1 msig implies 
they need a relatively large fraction of their stellar mass 
in older populations, and in all cases we find best fits 
with /old > 0.93. 

Chen & Marzke (2004) used a NICMOS-selected cata- 
log in the HUDF combined with photometric redshifts to 
select six candidates for massive, red (^775 — i^ieo ^ 2) 
early-type galaxies at z > 3. Interestingly, for only 2 
of their candidates where we have IRAC data does our 
analysis reach the conclusion that old stellar populations 
dominate the observed colors (NIC IDs 223, 625). We 
find for three of their candidates (NIC IDs 219, 1078, 
1237) that when the IRAC data are included the SED 
fits favor solutions where most of the light originates from 
dusty, young stellar populations. For the remaining can- 
didate of Chen & Marzke, we do not have IRAC data but 
our SED fit to the existing ACS -h NICMOS photometry 
does favor a solution where the galaxy is dominated by 
older stellar populations. 

For three of our six objects with IRAC photometry, 
z > 2.7, and Ud > 0.9 (NIC IDs 269, 989, and 1211), 
Chen & Marzke (2004) do not include these in their sam- 
ple of evolved galaxies. We derive a photometric redshift, 
z = 2.7, for NIC ID 269, slightly below the z > 2.8 selec- 
tion of Chen & Marzke, which may be the reason they 
excluded it. Chen & Marzke interpret object NIC ID 
1211 to be dominated by a dusty starburst. However, 
the HiQo to [3.6/im] and [4.5/im] colors imply that while 
its rest-frame UV light stems from dust-obscured young 
stars, most of the mass (>88% at 68% confidence) re- 
sides in evolved stellar populations (see figure 8; however, 
we note that the upper limit on the X-ray detection of 
this galaxy may imply that it hosts an AGN, see § 3.2). 
NIC ID 989 appears to have some ongoing star formation 
(with relatively blue ^435 to zg^o colors), but older stars 
appear to dominate the rest-frame optical and near-IR 
light (figure 7). We suspect this object was excluded 
from the selection of Chen & Marzke because of its rela- 
tively blue ACS colors. Regardless, we find that data at 
wavelengths greater than rest-frame 1 /im are required 
for interpreting the properties of the stellar populations 
of red galaxies. 

Parenthetically, we note that there are two objects in 
the HUDF with Jno - i^ieo > 2 mag (NIC IDs 1078, 
1211) and one object has Jno— ^leo = 1-95 mag (NIC ID 
1237). All three have IRAC photometry. These objects 
have colors comparable to the "unusual" red, Jno — ^i60- 
selected object in the HDF-N reported by Dickinson et 
al. (2000), which has Jno — Hiqq :^ 2.3 mag. Based on 
our fits to the ACS, NICMOS, and IRAC photometry of 
the HUDF objects, we conclude that two (NIC IDs 1078, 
1237) are dominated by dusty starbursts, while only one 
(NIC ID 1211) is dominated by a substantial population 
of old stars. Dickinson et al. (2000) are unable to distin- 
guish between these two possibilities (or, for that mat- 
ter, the possibility that their object is a z ?^ 12 Ly man- 
break- type galaxy). This is partly because the HDF- 
N object is undetected in the HST /WFFC2 data, nor 
did it have longer wavelength data from Spitzer /IRAC. 
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Two of the three objects here have HST/ ACS photom- 
etry > 1 mag fainter than the detection hmit of the 
WFPC2 data in the HDF-N. Therefore, we conclude 
that one requires very deep HST optical imaging and 
Spitzer /IRAC imaging to constrain the nature of these 
very red, Jno — i^ieo ^ 2 mag objects. 

None of our six candidates for galaxies with substan- 
tial population of evolved stellar stars is entirely devoid 
of recent star formation (i.e., all objects have /qm < 1-0). 
Indeed, in the entire HUDF there are no candidates for 
galaxies consistent with a purely passively-evolving stel- 
lar population formed at much higher redshift (and in 
the NICMOS HDF-N, there is only one possible candi- 
date, see discussion above; Dickinson et al. 2000). There- 
fore, such objects are extremely rare if they exist at this 
redshift. We conclude that at z ^ 2.5 even galaxies 
with a substantial amount of evolved stellar populations 
still maintain some ongoing star-formation. Papovich et 
al. (2005) argue based on galaxies' morphologies that at 
z > 2, recurrent star- formation dominates the rest-frame 
UV and blue light, which will maintain some level of ho- 
mogeneity in the galaxies' internal colors. Our interpre- 
tation here is consistent with that scenario. However, we 
note that it appears that while the young stellar popula- 
tions dominate the UV and blue light, in a small fraction 
of galaxies a substantial population of older stars domi- 
nates at rest-frame wavelengths > 1 fim. 

4.2. Galaxies with Evolved Stellar Populations at 
z > 1.5 

With the Jiio — i^ieo > 1 niag selection, we iden- 
tify massive galaxies with substantial stellar mass and 
with SEDs dominated by the light from passively evolv- 
ing stellar populations. Figure 9 shows the measured 
HiQo magnitudes as a function of redshift for galax- 
ies with /old > 0.9, for objects with IRAC photometry 
(right panel) and without IRAC photometry (left panel). 
Although the data should be sensitive to such objects 
down to our selection limit, i^ieo < 26 mag, only one 
such object has 25 < i^ieo < 26 mag — the rest have 
HiQo < 25 mag. Figure 10 shows the derived stellar 
masses as a function of galaxy redshift for all galaxies 
in the Jno — Hiqq > 1 mag sample. The figure shows 
the stellar masses derived for the whole sample using 
only the ACS through NICMOS data (left panel), and it 
shows the stellar masses derived including the IRAC data 
for the IRAC-detected subsample (right panel). Objects 
with more than 90% of their stellar mass in evolved stel- 
lar populations are circled in the figure. There exist 
objects whose light (and stellar mass) is dominated by 
passively evolving stellar populations with total stellar 
masses up to and exceeding >10^^ Mq over the entire 
redshift range. The existence of these galaxies is not 
entirely unexpected, given that such objects likely domi- 
nate the stellar mass density even at 2: ~ 2 (e.g., Rudnick 
et al. 2006; van Dokkum et al. 2006). 

However, there is an apparent dearth of objects whose 
light is dominated by passively evolving stellar popula- 
tions with lower stellar masses. For example, none of 
the 10 IRAC-detected galaxies with 2;phot > 1-6 and 
with more than 90% of their stellar masses evolved stel- 
lar populations are fainter than Hiqq > 25 mag (see 
figure 9), which is one magnitude above our detection 
limit. (We note there is one object with 2:phot — 1-5 



and H = 25.5 mag, and a derived total stellar mass 
4 X 10^ A^0, see figure 9.) Therefore, at these high red- 
shifts, the candidates for passively evolving galaxies are 
already massive. In terms of total stellar mass, the right- 
hand panel of figure 10 shows that all the IRAC-detected 
objects with /qm > 0.9 and z > 1.5 have stellar masses 
greater than ~2 x 10^^ A4q. If we instead consider the 
results for galaxies from their stellar masses derived with- 
out the IRAC data, there is only a single candidate as a 
low-mass galaxy whose mass is dominated by passively 
evolved stellar populations (and we consider the stellar 
masses derived in the old component without IRAC data 
to be less robust). 

The data should be sensitive to objects dominated by 
old stellar populations with stellar masses with ^^3 x 10^ 
to 10^° A^0(see figure 4 and figure 10). We note that this 
is not a bias related to the sources with IRAC detections. 
The IRAC data from GOODS (5.8 fim flux density of 
0.11 //Jy, 5(t) would detect a passively evolving stellar 
population formed at z=7 and M > 10^-^ Mq for all 
redshift s considered here. 

Any plausible mass function for galaxies at z ~ 2 
would predict a greater number of galaxies with lower 
stellar masses than at high stellar mass. For example, 
using the Fontana et al. (2006) empirical galaxy mass 
function, we would expect roughly double the number of 
galaxies with 0.5 — 2.0 x 10^^ Mq compared to those 
with >2 X 10^^ Mq. This strongly contrasts with the 
findings here for the red galaxies dominated by passively 
evolving stellar populations. Therefore, the data suggest 
that galaxies with more than 90% of their stellar mass in 
evolved stellar populations have a minimum mass thresh- 
old of roughly several times lO^^7W0atz~2. The lack 
of such objects with lower masses may be related to the 
processes governing stellar mass build-up in galaxies at 
these redshifts. 

4.3. Constraints on the Densities of Galaxies with 
Evolved Stellar Populations at 2.5 < z < 3.5 

With the full multiwavelength dataset, we have sep- 
arated our sample into those galaxies that are domi- 
nated by dust-enshrouded starbursts and those galaxies 
whose rest-frame optical and near-IR emission is consis- 
tent with being dominated by a substantial population 
of evolved stars (see, e.g., section 4.2, above). Of the 16 
Jiio — -f^ieo > 1 niag galaxies with robust IRAC photom- 
etry, 6 have best-fit models favoring a substantial pop- 
ulation of evolved stellar populations, which is 37.5% of 
the sample. Extending this to the full sample of 32 galax- 
ies with Jiio — HiQo > 1 mag, we estimate there should 
be 12 such objects in the area of the HUDF. Using these 
numbers, we place constraints on the number and mass 
density for these types of objects. If some of our 6 can- 
didates for evolved stellar populations turned out to be 
dusty starbursts, our numbers would be reduced. Due 
to the lower mass to light ratios of younger models, the 
fitted masses would also be reduced. Therefore the den- 
sities derived here should be regarded as upper limits. 

The HUDF data reach very deep limits in stellar mass 
for galaxies dominated by evolving stellar populations. 
For a galaxy whose stellar mass formed in situ at = 7, 
the HiQo = 26.0 mag detection limit corresponds to a 
limit in stellar mass ranging from 6 x 10^ M.q dit z = 2 
to 1 X 10^^ Mq d^t z = 4 (see figures 4 and 10). Although 
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Fig. 9. — -f/^ieo magnitudes vs. photometric redshifts for the subset of (Jno — ^16o)ab > 1-0 mag objects with best-fit /qm > 0.9, 
for objects without IRAC photometry (left panel) and with IRAC photometry (right panel). The dashed line shows the ifieoAB < 26 
magnitude cut. In the right-hand panel, cyan stars indicate those sources among the IRAC-detected subsample with /oid > 0.9 which have 
24 fim detections. Thick-circles denote objects with /oid > 0.9 and 2;phot > 2.5 



the (Jiio — HiQo) > 1 mag color selection should be, in 
principle, sensitive to old stellar populations at z > 2, 
in practice we find candidates for galaxies dominated 
by passively evolving stellar populations /om > 0.94 at 
z > 2.7. Likewise, at z > 3.5, the Balmer/4000A-breaks 
move beyond the F160W bandpass, making our color se- 
lection less efficient. Therefore, to set an upper limit on 
the number and mass densities of galaxies dominated by 
passively evolving stellar populations, we limit our red- 
shift range to 2.5 < z < 3.5. We note that requiring 
an IRAC detection does not limit our stellar mass selec- 
tion; the limiting stellar mass for the IRAC GOODS data 
sensitivity limit is always less than that for the NICMOS 
limit at these redshifts (see above). 

The volume in the area of the HUDF is 1.9 x 10"^ Mpc^ 
between 2.5 < z < 3.5. Restricting our analysis to only 
those galaxies with 2.5 < z < 3.5 we derive a number 
density of 3.1x10-'^ Mpc"^, or 6.3x10-"^ Mpc"^ includ- 
ing our 50% incompleteness. However, these six objects 
are not visible over the whole redshift range. In the num- 
ber density calculation, we compute the effective volume, 
Vef[ for each galaxy from z = 2.5 to Zq^^^ where Zeff 
is the lesser of z = 3.5 and the maximum redshift to 
which the galaxy could be detected. The number den- 
sity derived by summing T^il/vj^ over all i galaxies gives 
3.3l]^"5 X 10~^ Mpc~^, and after accounting for incom- 
pleteness, n = e.etg o X 10~^ Mpc"^. 

Assuming these objects will evolve to become present- 
day, passively evolving red galaxies, then this number 
density is significantly lower than what we would expect 
using local mass functions. Using the z ^ mass func- 
tion parameters listed in table 4 of iBell et al.l ( ^0 3), we 
expect the number density of red, early-type galaxies to 



be n{z = 0) = 3.5 x 10~^ Mpc~^ in the volume spanned 
between 2.5 < z < 3.5 in the HUDF (and again taking 
the redshift-dependent mass limit from figure [4j)^. This 
is larger than the observed number density (including in- 
completeness) by roughly a factor of 5. Therefore, while 
we would have expected 65 objects in the HUDF (i.e., 32 
after our 50% completeness) we have found only 6. 

This difference is substantially larger than the Pois- 
son noise (±a/6/6, i.e., 40%). It may be that these ob- 
jects suffer strong field-to-field clustering (possibly ex- 
pected for red galaxies, e.g., Daddi et al. 2003; Quadri 
et al. 2007). We estimate the number error from den- 
sity variations in the volume between 2.5 < z < 3.5 in 
the HUDF as 0.3 x ag. While we have no measurement 
for (78 for these gal axies, based on theore tical models we 
expect as ~ 1.0 ( Adelberger et al1l2005f ), implying that 
the error from density variations is 30%. Summing these 
errors in quadrature brings our total error estimate to 
^50% for the expected number of these galaxies in the 
HUDF. Even so, the difference between the 32 expected 
galaxies and the 6 we find is substantially larger than 
the combined errors from counting statistics and cosmic 
variance. Therefore, by number, and recalling that these 
candidates only define upper limits, the objects in the 
HUDF which are dominated by evolved stellar popula- 
tions make up less than one-third of present-day red, 

^ Note that fBell et al."2003') assumed a formation redshift, Zf = 
4, in their mass functions. This is only marginally lower than 
Zf = 7 used in this analysis. However, the change in lookback time 
from 2; ~ to 4 relative to that from 2; ~ to 7 corresponds to a 
change in the mass-to-light ratio of M/Lb = 5.8 to 6.2. Applying 
our higher formation redshift to Bell et al. would increase their 
stellar masses by ~ 5%, which is negligible compared to other 
uncertainties. 
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Fig. 10. — Best-fit model masses and photometric redshifts for all HUDF objects with (Jno — ^^16o)ab > 1-0 mag (left panel) and the 
subset of these with unblended IRAC photometry (right panel). The dashed line shows the limiting stellar mass for galaxies with a stellar 
population formed in a burst at z = 7 with observed magnitude ifieo = 26.0. In both panels the open circles indicate objects dominated 
by old stellar populations (/old > 0-9) • In the right-hand panel, cyan stars indicate those sources among the IRAC-detected subsample 
with 24 /im detections. Thick-circles denote those objects with /qM > 0.9 and 2;phot > 2.5, which are the discussion of § 4.3. 



early- type galaxies. Thus, at least two- thirds of present- 
day red, early- type galaxies must assemble and/or evolve 
into their current configuration at z < 2.5. 

Much of the stellar mass density in galaxies 2.5 < z < 
3.5 could reside in red-selected galaxies, whose light is 
dominated by older, passively evolving stellar popula- 
tions. For the Juq—Hiqo > 1 mag galaxies in the HUDF, 
we derive an upper limit on the stellar mass density of 
log(p7[M0Mpc"^]) = 7.4 ± 0.3, including the effect of 
our redshift-evolving mass limit, where our errors are de- 
rived from a Monte-Carlo simulation of the data. When 
we include our 50% incompleteness, our total stellar mass 
density becomes \og{p* /[MqMpc~^]) = 7.7 ± 0.3. Pois- 
son sampling of the 6 objects gives a factor of two error 
and therefore is the dominant source of uncertainty in 
the mass density estimate. We consider this measure- 
ment to be an upper limit on the mass density of distant 
evolved galaxies for two reasons: a) it is possible that 
some of our evolved galaxy candidates are in fact dusty 
starbursts, and b) because our two-component model fit 
incorporates an a priori upper limit on the mass-to-light 
ratio in the form of a maximally old SSP component (see 

iRudnick et al.l (|2003l ) derive a stellar mass density at 
z = 2.S0toig of log P* = ^'^-o.i and at z = 2.0ll^-j? 
of log = 7.5^0 1- Because our mass density provides 
an upper limit on the mass limit, we conclude that our 
derived mass density is consistent those of Rudnick et 
al. In this case, it may be that passively stellar popula- 
tions dominate the stellar mass density. This is consis- 
tent with the analysis of Dickinson et al. (2003) who used 
two-component star-formation history models to set an 



upper limit on the stellar mass density from passively 
evolving stellar populations hidden beneath the glare of 
younger stars in galaxies at z < 3. Dickinson et al. find 
that at z ~ 2.5 — 3 an upper limit on the stellar mass den- 
sity of log(p*/[A^0Mpc"^]) = 7.89t^;?^. If this reflects 
reality, then based on our analysis, almost all of this 
mass density coul d reside in faint, red -selected galaxies 
at these redshifts. Ivan Dokkum et al] p006) found that 
at z ~ 2 red-selected galaxies with masses greater than 
10^^ Mq dominate the mass density. Our analysis of 
the HUDF galaxies suggests that at lower stellar masses 
red-selected galaxies still contribute significantly to the 
total mass budget. 

Lastly, we stress that with only rest-frame UV and 
optical data, it is not possible to discriminate between 
those galaxies dominated by dusty starbursts and older 
stellar populations (e.g., Moustakas et al. 2004; Small et 
al. 2002). Although the IRAC data does not resolve this 
degeneracy for all galaxies, we show that a sub-set of ob- 
jects do have consistent best-fit results and we therefore 
use these objects in our analysis. 
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